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ABSTRACT 

We study the evolution of a protoplanct of a few earth masses embedded in 
a protoplanetary disc. If we assume that the atmosphere of the protoplanet, 
i.e. the volume of gas in hydrostatic equilibrium bound to the core, has a 
surface radius smaller than the Roche lobe radius, we show that it expands as 
it accretes both planetesimals and gas at a fixed rate from the nebula until it 
fills in the Roche lobe. The evolution occurs on a timescale shorter than the 
formation or migration timescales. Therefore, we conclude that protoplanets of 
a few earth masses have an atmosphere that extends to the Roche lobe surface, 
where it matches to the nebula. This is true even when the Bondi radius is 
smaller than the Roche lobe radius. This is in contrast to the commonly used 
models in which the static atmosphere extends up to the Bondi radius and is 
surrounded by a cold accretion fiow. As a result, any calculation of the tidal 
torque exerted by the disc onto the protoplanet should exclude the material 
present in the Roche lobe, since it is bound to the protoplanet. 

Key words: planetary systems — planetary systems: formation — planetary 
systems: protoplanetary discs — planets and satellites: general 



1 INTRODUCTION 

The recent release of Kepler data reporting more than 1200 planet candidates transiting their 
host star confirms that the core accretion scenario for forming planets is ubiquitous. Indeed, 
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most of these objects are less than half the radius of Jupiter, which indicates that they are 
at most Neptune-size (Borucki et al. 2011). Such low-mass objects are probably not the 
result of gravitational instabilities, but are more likely formed through the accumulation of 
a metal rich core followed by the capture of a more or less massive envelope of gas (Lissauer 
1993 and references therein). 

In this model, the core is beheved to be formed through the solid-body accretion of 
kilometre-sized planetesimals. For typical disc parameters, when the core reaches about 0.1 
earth mass, it starts binding the gas of the nebula in which it is embedded and a gaseous 
atmosphere forms around it. As long as the core is less massive than the so-called criti- 
cal mass, the energy radiated from the atmosphere is compensated for by the gravitational 
energy that the planetesimals entering the atmosphere release when they collide with the 
surface of the core. During this stage of evolution, the atmosphere is therefore in quasi-static 
and thermal equilibrium and grows slowly in mass along with the core (Perri & Cameron 
1974, Mizuno 1980). By the time the core reaches the critical mass, the atmosphere has 
become too massive to be supported at equilibrium by the energy released by the planetesi- 
mals. At that point, it starts contracting and the subsequent runaway accretion of gas leads 
to the formation of giant planets (Bodenheimer & Pollack 1986, Pollack et al. 1996). 

At the same time as they form, protoplanets migrate through the nebula as a result of 
tidal interaction with the surrounding gas (Goldreich and Tremaine 1979, 1980, Lin and 
Papaloizou 1979, 1993 and references therein, Papaloizou & Lin 1984, Ward 1986, 1997). 
Calculations performed in isothermal discs show that cores of several earth masses migrate 
inward on a relatively short timescale, shorter than the planet formation timescale (Ward 
1997, Tanaka et al. 2002, Bate et al. 2003). However, more recent calculations including 
detailed energy balance in non isothermal discs show that migration can be slower or even 
outward (Paardekooper & Mellema 2006, Kley & Crida 2008, Paardekooper & Mellema 2008, 
Baruteau & Masset 2008, Paardekooper & Papaloizou 2008). Usually, the torque exerted by 
the disc on the protoplanet is calculated by excluding the gas comprised in its Roche lobe. 
However, it has been suggested by D'Angelo et al. (2003, see also Crida et al. 2009) that the 
gas present in the Roche lobe but not bound to the protoplanet may contribute significantly 
to the total torque and slow down migration. This happens when the atmosphere of the 
protoplanet, i.e. the volume of gas at hydrostatic equilibrium that is bound to the core, does 
not fill in its Roche lobe, and is surrounded by a cold (Bondi type) accretion flow. 

Models of atmospheres of protoplanets that do not fill in the Roche lobe have been con- 
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sidered by different authors. Papaloizou & Nelson (2005) have constructed models assuming 
that the only source of energy for the atmosphere is that produced by the gravitational 
contraction of the gas. This corresponds to the runaway gas accretion phase, i.e. after the 
core has reached the critical mass, as the luminosity due to the accretion of planetesimals 
is then negligible. The radius of the protoplanet atmosphere may become smaller than the 
Roche lobe radius in this phase if the protostellar disc cannot supply gas to the atmosphere 
rapidly enough. Such models have also been considered by Mordasini et al. (2011). Lissauer 
et al. (2009, see also Movshovitz et al. 2010) have calculated models of atmospheres prior 
to the runaway gas accretion phase by assuming that the surface radius is either the Bondi 
radius or a fraction of the Roche lobe radius depending on whether the Bondi radius is 
smaller or larger than the Roche lobe radius. This is based on the results of three dimen- 
sional numerical simulations indicating that only gas within about one quarter of the Roche 
lobe remains bound to the protoplanet. The gas accretion rate onto the atmosphere is then 
taken to be the value required to match the outer radius to the desired value. To calculate 
the volume of gas which is bound to the protoplanet from the numerical simulations, they 
identify the trajectories of tracer particles that are trapped inside the gravitational potential 
of the protoplanet. We comment that this approach may severely underestimate the amount 
of gas that can be accreted as in reality particles can collide with the protoplanet and be- 
come bound if energy is dissipated into shocks, as in the process of star formation. If the 
protoplanet fills in a significant part of its Roche lobe, collisions may become frequent. 

Three dimensional numerical simulations of protoplanets of a few earth masses in proto- 
planetary discs usually start with a core surrounded by a static atmosphere that extends at 
most up to the Bondi radius. The Keplerian flow located within and around the Roche lobe, 
and in which the protoplanet is embedded, is then perturbed and material from the disc 
rains down on the protoplanet from above and below, because the gravitational potential 
from the protoplanet is not balanced (D'Angelo et al. 2003, Paardekooper & Mellema 2008, 
Machida et al. 2010). This results in the protoplanet being surrounded by a cold isother- 
mal accretion flow. In the present paper, we show that, if a core below the critical mass is 
surrounded by an atmosphere that does not fill in entirely the Roche lobe and is embedded 
in an accretion flow, it is going to evolve into a static model expanded to the Roche lobe 
radius on a very short timescale, even when the Bondi radius is smaller than the Roche lobe 
radius. 

In the calculations presented here, we include the accretion of planetesimals onto the core. 
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as we focus on the stage of the evolution before runaway gas accretion starts. In section |2], 
we calculate the mass of the atmosphere that can be supported at hydrostatic and thermal 
equilibrium around a core of a given mass assuming the surface radius of the atmosphere 
to be fixed. We consider a range of surface radii from a few core radii up to the Roche lobe 
radius. In section [3l we consider the evolution of a core that starts accreting gas at a given 
rate from a radius much smaller than its Roche lobe radius. In contrast to what is usually 
assumed, we suppose here that the gas accretion rate onto the protoplanet is fixed, as it 
is likely to be constrained by the disc's physical parameters (e.g., density, temperature). 
We calculate a sequence of quasi-static equilibria along which the mass of the protoplanet 
increases. We show that the atmosphere has to expand very rapidly to accommodate the 
gas that is being accreted, even if the accretion rate is relatively small. We argue that, even 
when the Bondi radius is smaller than the Roche lobe radius, the protoplanet will eventually 
fill in its Roche lobe. Finally, in section HJ we summarise and discuss our results. 



2 DEPENDENCE OF THE ATMOSPHERE MASS ON THE 
ATMOSPHERE SURFACE RADIUS 

We consider an atmosphere at hydrostatic and thermal equilibrium and calculate its mass 
for a fixed core mass, accretion rate of planetesimals and surface radius of the atmosphere. 

2.1 Structure of the protoplanet atmosphere 

The equations governing the structure of the protoplanet atmosphere at hydrostatic and 
thermal equilibrium are the same as the equations of stellar structure, with the nuclear en- 
ergy rate being replaced by the rate at which planetesimals that enter the planet atmosphere 
and collide with the core surface release their energy. These equations have been presented 
in Papaloizou & Terquem (1999) and we recall them below. 

We assume that the protoplanet is spherically symmetric and nonrotating. We denote r 
the radius in spherical coordinates in a frame with origin at the centre of the protoplanet. 
The equation of hydrostatic equilibrium is: 
dP 

= -'P- 

Here, P is the pressure, g = GM{r)/r'^ is the acceleration due to gravity, with M(r) being 
the mass interior to radius r (this includes the core mass if r is larger than the core radius) 
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and G is the gravitational constant. The relation between M(r) and the mass density per 
unit volume p is: 

— = 47rrV- (2) 

dr 

We adopt the equation of state for a hydrogen and helium mixture given by Chabrier 
et al. (1992) for mass fractions of hydrogen and helium of 0.7 and 0.28, respectively. The 
standard equation of radiative transport gives the luminosity Lrad that is transported by 
radiation through the atmosphere as a function of the temperature gradient dT/ dr: 

dT _ -3k,P Lrad /ON 

dr ~ WaT^Anr^' ^ ' 

where k is the opacity, which in general depends on both p and T, and a is the Stefan- 
Boltzmann constant. 

As a large part of the atmosphere interior is unstable to convection, only a fraction of the 
total luminosity is transported by radiation. In the models presented here, we assume that 
the only energy source comes from the planetesimals that are accreted by the protoplanet 
and release their gravitational energy as they collide with the surface of the core. They 
generate a total core luminosity L^ove given by: 

-•^core — ) V / 

where Mf. and Tc are, respectively, the mass and the radius of the core, and Mc is the 

planctcsimal accretion rate. 

The radiative and adiabatic temperature gradients, Vrad and Vad, are given by: 

_ / d\nT \ _ 3/tLcoreP f s 

~ \dhiP j^^ ~ 647raGMT^' 

and 

with the subscript s denoting evaluation at constant entropy. 

In the regions where Vrad < Vad, there is stability to convection and therefore all the energy 
is transported by radiation, i.e. Lrad — -t'core- In the regions where Vrad > Vad, there is 
instability to convection and therefore part of the energy is transported by convection, i.e. 
Lcore — Ly-ad + Lconv, where Lconv is the luminosity associated with convection. Using the 
mixing length theory (Cox & Giuli 1968), wc obtain: 



-^conv — '^f CpA^Y 



'dT\ _ /dT' 
dr I \ dr 



- 3/2 


1 






^P9 





(7) 



6 Terquem and Heinemann 

where Ami = \cimiP/{dP/df^)\ is the mixing length, ami being a constant of order unity, 
(dT/dr)^ = VadT (dlnP/dr), and the subscript P denotes evaluation at constant pressure. 
The different thermodynamic parameters needed in the above equation are given by Chabrier 
et al. (1992), and we fix ami = 1- 

2.2 Boundary conditions 

We have to solve the above equations for the three variables P, M and T as a function of 
r. Accordingly, we need three boundary conditions. 

Following Bodenheimer & Pollack (1986) and Pollack et al. (1996), we take for the mass 
density of the core pc = 3.2 g cm~^. The inner boundary of the atmosphere is equal to the 
core radius Vc, given by: 



(8) 



The first boundary condition is that M(rc) = Mc- 

We denote ratm the outer boundary of the atmosphere, and we assume this radius to be 
fixed. In the calculations presented below, we will take ratm to be either a few core radii or 
the Roche lobe radius r^, which is assumed to be the same as the Hill radius, and is given 
by: 

where a is the orbital radius of the protoplanet and is the mass of the central star. 
Throughout the paper, we take = 1 Mq. 

When Tatm < ^L, we assume that the space above the atmosphere and in the Roche 
lobe of the protoplanet is filled with a cold accretion flow. Then the temperature Tg and 
the pressure Ps at the surface r = ratm can be approximated by the standard photospheric 
conditions (e.g., Kippenhahn & Weigert 1990, Bodenheimer et al. 2000): 



L 



1/4 



where a is the Stefan-Boltzmann constant, and: 



r 



where Kg = n {ps, Tg) and ps is the mass density at r = ratm- Similar boundary conditions have 
been used by Palla & Stabler (1992) for computing the evolution of a protostar accreting 
from a disc. The validity of these boundary conditions will be discussed in section 13. 1[ 
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Note that in any case the structure of the atmosphere does not depend significantly on Tg 
(Papaloizou & Terquem 1999). 

We also present below some calculations with ratm = 'Tl, which are done using the 
boundary conditions presented in Papaloizou & Terquem (1999). 



2.3 Calculations 

For a given core mass M^, planetesimal accretion rate and atmosphere radius ratm? we 
solve the equations ([1]), ([2]) and ([3]) with the boundary conditions described above to get the 
structure of the atmosphere and its mass Matm- We start the integration at r = ratm with 
T = Ts, P = Ps and an initial guess for M(ratm) = + ^atm- In principle, to calculate 
Ps, we need to know ps, which we derive from the pressure. However, for the values of Tg 
we obtain, which are smaller than 1000 K, the opacity Kg does not depend on pg. We can 
therefore self-consistently ignore pg when calculating Kg. The integration is carried through 
down to r = r^. In general, the condition that M(rc) = will not be satisfied. An iteration 
procedure is then used to adjust the value of M(ratm) until M(rc) = to a specified 
accuracy. When the core mass is larger than the critical mass Mcnt, no solution can be 
found, i.e. no atmosphere at equihbrium can be supported at the surface of a core which 
mass exceeds Merit. 

The opacity is taken from Bell & Lin (1994). This has contributions from dust grains, 
molecules, atoms and ions. It is written in the form k = nip°'T^ where Kj, a and h vary with 
temperature. 

We take Mc to be in the range 10~*-10~^ yr~^, as this gives timescales for building 
up a 10 core of 10^-10^ years. These values are consistent with the calculation of M^. by 
Pollack et al. (1996) and, more recently, by Movshovitz et al. (2010). 

Figure [1] shows the mass of the atmosphere versus that of the core for = 10~^, 
10~^, 10~^ and 10~^ yr~^, and for an atmosphere radius ratm equal to Sr^, lOr^ or 
r^. When ratm = 5rc, the photospheric temperature and pressure are Tg = 555.3 K and 
Pg = 3.73 erg cm'^ for = 10"^ M® yr-\ and Tg = 98.7 K and Pg = 20.42 erg cm^^ 
Mc = 10~® yr^^- The curves corresponding to ratm = are calculated assuming the 
temperature and pressure in the midplane of the nebula to be 33 K and 0.025 erg cm~^, 
respectively, which correspond to a standard steady state disc model with a = 10~^ and gas 
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accretion rate M = 10~^ Mq yr~^, assuming the protoplanet is at 5 AU from the central 
star (see Papaloizou & Terquem 1999). 

For Tatm = lOrc, the critical core mass is 36, 24, 16 and 10 for Mc = 10~^, 10~^, 
10^'' and 10^*^ yr^"*^, respectively. For ratm = Sr^, the critical core mass is 55, 37, 25 and 
16 for these values of Mc. We note that, if the core is at 5 AU from the central star, 
'"atm = lOrc and ratm = represent about 1% and 0.7%, respectively, of the Roche lobe 
radius of the protoplanet when its core reaches the critical mass. 

For comparison, in a standard steady state disc model with a = 10~^ and gas accre- 
tion rate M = 10~^ Mq yr~^, at 5 AU from the central star, the critical core mass of a 
protoplanet with ratm = is 24, 15, 9 and 5 for = 10"^ 10-^ 10"^ and 10"^ 
Mq yr~^, respectively. These values are almost independent of the conditions in the nebula, 
and therefore do not change when M and a are varied. 

We see that the critical core mass increases when ratm decreases. This is because the gas 
in the atmosphere has less (negative) gravitational energy when ratm is smaller, i.e. it is more 
bound, and therefore, for a fixed core accretion luminosity, it can be supported at equilibrium 
up to larger core masses. Also, the temperatures at the bottom of the atmosphere tend to 
be lower for smaller ratm, so that more mass can settle near the surface of the core, which 
also helps to increase the critical core mass. Finally, because the surface of the atmosphere 
is smaller, radiative losses are smaller which helps supporting the atmosphere at equilibrium 
up to larger core masses. However, for a given core mass and core accretion luminosity, the 
mass of the atmosphere decreases as ratm decreases. 

In all the calculations displayed on figure [H, ri is smaller than the Bondi radius tb and 
represents the largest radius the static atmosphere of the protoplanet can expand to. If 
the midplane temperature of the protoplanetary disc were larger, r^ would become smaller 
than r^. This is the case, for instance, for protoplanets less massive than 16 Mq when the 
midplane temperature is 140 K, which corresponds to a a disc model with a = 10~^ and gas 
accretion rate M = 10"'^ Mq yr~^, at 5 AU from the central star. 

3 EVOLUTION OF A PROTOPLANET ACCRETING FROM A 
PROTOPLANETARY DISC 

We are now going to consider the evolution of a protoplanet which has an atmosphere that 
does not fill in its Roche lobe. 
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3.1 Accretion onto a core 

Numerical simulations of low mass protoplanets (~ a few Mq) embedded in protoplanetary 
discs and which do not fill entirely their Roche lobe show that they significantly perturb the 
Keplerian motion of the fluid in the vicinity of and within their Roche lobe. High resolution 
two dimensional computations by D'Angelo et al. (2002) indicate that circumplanetary discs 
may form within the Roche lobe of cores more massive than about 5 M^. Subsequent three 
dimensional simulations by the same authors (DAngelo et al. 2003; see also Bate et al. 2003), 
using an isothermal equation of state, show that the circumplanetary disc features are much 
less marked around low mass cores due to the vertical motion of the fluid around them. Gas 
entering the Roche lobe from above and below the protoplanet falls down directly onto it 
instead of joining a circumplanetary disc in the midplane. However, in these calculations, 
the resolution is limited to be several percent of the Hill radius of the protoplanet, which is 
between one and two orders of magnitude larger than the radius of the core itself. Therefore, 
from these simulations, the formation of a circumplanetary disc around low mass cores cannot 
be ruled out. Simulations by Machida et al. (2008) and Machida (2009) of circumplanetary 
disc formation around protoplanets show indeed that these discs tend to be rather compact, 
with an outer radius smaller than about 50 core radii. 

Further three dimensional numerical simulations by Paardekooper & Mellema (2008), 
including radiative transfer, show significant differences compared to the isothermal case. In 
particular, the accretion rate onto the protoplanet is much lower due to the compression of 
gas within the planetary atmosphere (which, in these calculations, is the unresolved region 
around the core) . The authors comment that the formation of a hot bubble around the pro- 
toplanet limits the gas flow towards its surface. In these simulations, the resolution is again 
on the order of 1% of the Hill radius, so that the potential formation of a circumplanetary 
disc close to the core could not be seen or dismissed. 

We suppose here that the protoplanet accretes through either a cold Bondi type accretion 
flow or a circumplanetary disc. The surface of the protoplanet atmosphere is at a radius Tatm 
which is the transition from the (spherical or disc-like) accretion flow to a hydrostatic 
model. The value of the gas accretion rate Mgas onto the protoplanet is very uncertain. The 

simulations by DAngelo et al. (2003) indicate an infall rate onto a 5 core on the order 
of 10~^ yr~^, but computations incorporating radiative transport give a value which is 
an order of magnitude smaller (Paardekooper & Mellema 2008). Computations with higher 
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resolution would give even smaller values as compression of the gas near the core would be 
larger on smaller scales. 

When the gas accreted encounters the surface of the protoplanet atmosphere, it releases 
the luminosity Lace = G {Mc + Matm) ^gas /''"atm- We assume that essentially all this lumi- 
nosity is converted into heat in and radiated away from the accretion shocks that form at the 
surface of the protoplanet (we neglect the energy that is converted into rotational energy of 
the protoplanet). The gas that is accreted though has to bring some entropy deep inside the 
atmosphere, since the temperature and pressure rise there as the atmosphere mass increases, 
and it does so by inducing some gravitational contraction of the atmosphere. This results 
in an increase of the internal energy of the atmosphere which in turn drives the expansion 
of its surface radius. However, since the mass of gas which is accreted is brought to the 
surface of the atmosphere, the resulting increase in the protoplanet luminosity is not very 
important (see discussion in section 7.2 of Papaloizou & Nelson 2005). We will neglect it 
here and will therefore consider that the accretion of gas does not modify the luminosity of 
the protoplanet, which is always given by Lcore- 

In section 12.21 we have assumed that the temperature and pressure at the surface of 
the static atmosphere were given by the photospheric values. This is a valid approximation 
as long as the temperature and pressure of the accretion flow are low compared to the 
photospheric values for the planet models. Note that if some of the accretion luminosity 
were not radiated away, which may well be the case, the surface temperature Tg would be 
higher than the photospheric value. However, we remark that the structure of the atmosphere 
is not very sensitive to for the masses considered here and, in any case, an increase of 
Ts would make the atmosphere more extended and therefore would strengthen the point we 
make in this paper. As far as the surface pressure is concerned, the only way the infalling 
gas contributes to it is through its ram pressure Pram, because its motion is supersonic 
(e.g., Paardekooper & Mellema 2008). When the accretion is spherically symmetric, the ram 
pressure is given by: 

where we have assumed that the gas is accelerated to free-fall velocity by the gravitational 
potential of the protoplanet. In the calculations discussed in section 13.21 below, we have 
checked that Pram is always smaller (and in most cases much smaller) than the photospheric 
pressure Pg. 



2G(M, + Matm^T'/' 



(12) 
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How exactly the gas is accreted by the protoplanet is not known, but we assume that once 
the gas is accreted it is redistributed in a spherically symmetric way at the surface of the 
atmosphere. As we mentionned above, the gas that is accreted induces some gravitational 
contraction of the atmosphere. However, as the accretion rate is low, we can consider that the 
atmosphere stays at quasi hydrostatic and thermal equilibrium as it accretes. The timescale 
for establishing hydrostatic equilibrium is indeed very short. If the temperature everywhere 
in the atmosphere were equal to that in the midplane of the protoplanetary disc (this gives 
a lower limit on the real value), i.e. between 30 and 140 K for a typical disc at a distance 
of 5 AU from the central star, then the sound crossing time through a distance of 10^^ cm, 
which is a typical value of the Hill radius of a few earth mass protoplanet, would be on the 
order of a year. The amount of gas accreted on such a timescale is a tiny fraction of the 
mass of the atmosphere, and so it can be assumed that it is brought to equilibrium almost 
inst ant aneously. 

We now consider the evolution of the protoplanet as it accretes gas at constant luminosity, 
assuming quasi hydrostatic and thermal equilibrium. 

3.2 Evolution of a 5 protoplanet 

For illustrative purposes, we consider a 5 core accreting from a protoplanetary disc. 
We calculate a sequence of quasi-static equilibria along which the mass of the protoplanet 
increases. As the gas accretion rate is fixed, constrained by the physics in the disc, the surface 
radius ratm of the atmopshere has to adjust to accomodate the mass increase. If the mass 
of the protoplanet atmosphere, with initially ratm < fi, increases rapidly relative to that 
of its core, then, for a constant luminosity, wee see from figure [H that ratm has to increase, 
i.e. the atmosphere of the protoplanet expands. In that case, the evolution is along almost 
vertical lines going from the lower curve (ratm = Sr^ on these plots) to the upper curve, 
which corresponds to ratm = ^l- 

For all reasonable values of Mgas and M^, we have found that a 5 core always 
undergoes such an evolution and that its atmosphere has to fill in its Roche lobe after a 
rather short timescale. To illustrate this point, we now develop in detail some particular 
cases, starting with a core that has ratm = Sr^. 

If Mc = 10"^ yr"\ such a value of ratm corresponds to Matm = 3.3 x 10"^ M^. We 
assume that the atmosphere accretes at a rate Mgas = 10~^ yr~^ After 10^ years, the 
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core mass has hardly changed, but Matm has grown significantly, being about 1.3 x 10~^ M®, 
which corresponds to ratm — 30rc. After 4 x 10^ years, the atmosphere mass has grown to 
4.3 X 10~^ Mq, which correspond to ratm = rL- 

Let us now consider the case Mc — 10~^ yr~^, which gives Matm = 9.2 x 10~^ 
for ratm = 5rc, and Mgas = 5 x 10""^ yr~^ After 5 x 10^ years, Mc = 5.05 and 
Matm = 3.4 X 10~^ M®, which corresponds to ratm = 25rc. After 2.5 x 10^ years, = 5.25 
and Matm = 1-3 x 10^^ M®, which corresponds to ratm = tl- 

Only when or Mga^ get very small may we have models with ratm < f^L- This is 
the case, for instance, when Mc = 10~^ yr~^ and Mgas < 10~^ yr~^, or when 
Me = 10-^ Me yr-^ and Mgas < 5 x 10"^ yr-^ 

Therefore, the static atmosphere of a protoplanet of a few earth masses has to fill in 
the Roche lobe as gas is being accreted. There is no Bondi type accretion flow around the 
protoplanet. The rate at which gas penetrates from the nebula through the Roche lobe 
surface in the atmosphere depends on the rate at which the core grows. As the core grows in 
mass from the accretion of planetesimals, gas from the nebula is allowed to penetrate in its 
atmosphere to keep it at equilibrium. Unless the core grows very fast or the protoplanet is 
in a very low density part of the nebula, there will always be enough gas coming through the 
Roche lobe surface to supply the atmosphere so that it always fills in the Roche lobe. Only 
when the core reaches the critical core mass may the atmosphere contracts so that ratm < t^l- 
That happens if the nebula cannot deliver gas fast enough to the collapsing atmosphere. We 
therefore conclude that a protoplanet with a core mass below the critical mass always fills in 
the Roche lobe. 

Note that the point we are making in this paper would still be valid if the luminosity of 
the protoplanet were not generated by the accretion of planetesimals but by the contraction 
of the atmosphere, provided the Kelvin-Helmholtz timescale were long enough so that the 
atmosphere were at quasi-hydrostatic and thermal equilibrium over the accretion timescale. 

We have assumed here that the atmosphere can freely expand up to the Roche lobe 
radius. We now give an argument to show that this is always the case, even when the Bondi 
radius is smaller than the Roche lobe radius. 
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3.3 Roche lobe radius vs. Bondi radius 

It is usually assumed that a protoplanet embedded in a protoplanetary disc cannot expand 
beyond the Bondi radius defined as: 

= 2 

where c is the sound speed in the protoplanetary disc (e.g., Bodenheimer et al. 2000). The 
argument is that the gas located at a distance from the protoplanet larger than tb has a 
thermal energy larger than the gravitational energy that would bind it to the protoplanet, 
and therefore cannot be accreted by it. However, if a molecule of gas located in the pro- 
toplanetary disc beyond is accelerated toward the protoplanet and collides with it, it 
may become bound (and therefore be accreted) like any other molecule coming from within 
the Roche lobe as described in section 13.11 Accretion of this molecule only requires that 
it loses at least part of its (kinetic plus thermal) energy into shocks during the collision. 
If the molecule hits and settles into a circumplanetary disc before falling onto the planet, 
accretion will happen as energy is radiated away from the disc. When the molecule becomes 
bound to the protoplanet, the pressure gradient in the atmosphere adjusts itself to balance 
the gravitational attraction of the protoplanet, i.e. hydrostatic equilibrium is maintained. 
At fixed luminosity, the protoplanet has to expand to accommodate this extra mass, and 
its final surface radius does not have to be limited by r^. For a given mass M-^tm and sur- 
face temperature of the atmosphere, an equilibrium atmosphere can be constructed with 

'^atm > rB- 

This is exactly similar to the process of star formation. In a molecular cloud that collapses 
onto a star, when particles hit the surface of the star they have a positive energy. This 
is because they have some thermal energy in addition to the kinetic energy due to the 
acceleration by the gravitational potential of the forming star. Accretion is possible because 
at least part of this energy is dissipated into shocks. The particle that is accreted may or 
may not bring some entropy into the stellar envelope, depending on whether all or only some 
of the energy is radiated away, and the star readjusts to equilibrium (Palla & Stabler 1991, 
1992, Hosokawa et al. 2010). 

We conclude that the protoplanet can always expand to fill in its Roche lobe, even when 
the Bondi radius is smaller than the Roche lobe radius. 



(13) 
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4 DISCUSSION AND CONCLUSION 

We have shown above that, if a 5 earth mass protoplanet has a static atmosphere that does 
not fill in the Roche lobe but is embedded in a cold accretion fiow, it has to expand as it 
accretes both planetesimals and gas from the nebula until it does fill in the Roche lobe. The 
evolution occurs on a timescale shorter than the formation or migration timescales. The 
value of 5 earth masses was chosen to illustrate the point, but the conclusion of course holds 
for any smaller protoplanets. Therefore, we conclude that protoplanets of a few earth masses 
have an atmosphere with a static structure all the way up to the Roche lobe surface, where 
it matches to the nebula, and that remains the case until the core reaches the critical mass 
or the gas in the nebula gets depleted. 

This result holds whether or not the Bondi radius is smaller than the Roche lobe radius. 
Therefore the Bondi radius does not determine the extent of the static structure around the 
core, although it gives the scale on which the density increases in the atmosphere. For radii 
smaller than the Bondi radius, the structure of the atmosphere is almost independent of the 
conditions in the nebula, whereas beyond the Bondi radius the density does not vary much 
and the structure is determined by matching to the pressure and temperature in the nebula. 

Numerical simulations which show that low mass protoplanets affect the flow in the 
surrounding protoplanctary disc inside their Roche lobe start with artificial initial conditions. 
The flow in the Roche lobe of the protoplanet is not at equilibrium to start with, as the 
gravitational attraction of the core is not balanced. Therefore, as seen in the simulations, 
material from the disc rains down on the protoplanet (D'Angelo et al. 2003; Paardekooper 
& Mellema 2008). When radiative effects are taken into account, this infall is slowed down 
as the gas is being compressed around the protoplanet (Paardekooper & Mellema 2008, 
Ayliffe & Bate 2009). When the luminosity generated by this compression is high enough, 
an almost spherical atmosphere forms around low-mass protoplanets, ultimately filling in 
their Roche lobe (Ayliffe & Bate 2009). This is in agreement with the results presented in 
this paper, where the luminosity is not being produced by the gas compression but by the 
accretion of planetesimals. The results reported in the present paper indicate that numerical 
simulations of low-mass cores in discs, when they cannot result in the build-up of a static 
atmosphere around the core (like in the isothermal case), should have as initial conditions a 
core surrounded by a static atmosphere extending all the way up to the Roche lobe radius. 
In such a case, there will be no cold accretion flow around the protoplanet. 
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It has been argued that the orbital migration of low mass cores depends very significantly 
on the torque exerted by the gas within the Roche lobe of the protoplanet (D'Angelo et al. 
2003, Crida et al. 2009). The results presented above indicate that this is probably not the 
the protoplanet itself fills in its Roche lobe. 
In this paper, we have considered the evolution of a protoplanet which accretes planetes- 
imals from the protoplanetary disc. Such an accretion may be reduced significantly before 
the gas in the nebula has disappeared, with the consequence that the core reaches the critical 
mass and the atmosphere has to contract to provide the energy radiated at its surface. In 
that case, its evolution is that calculated by Papaloizou & Nelson (2005). 

Finally, we comment that we have not taken into account here the rotation of the proto- 
planet onto itself. Although this would not affect the main conclusion of our paper, we point 
out that there may be some boundary layer at the surface of the protoplanet's atmosphere 
due to a finite relative velocity (shear) between the atmopshere and the surrounding disk 
material. The size of the boundary layer would be determined by the disc's viscosity. 

We thank S. Balbus for discussions and J. Papaloizou for very useful comments which 
have greatly improved the original version of this paper. C.T. thanks Princeton Univer- 
sity, where the research reported in this paper was done, for hospitality and support. T.H. 
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Figure 1. Atmosphere mass, Matm, in units of Mg, vs. core mass, Mc, in units of Mg, for different planetesimal accretion 
rates Mc. The different frames correspond to Mc = 10~* (top left), 10^^ (bottom left), 10~^ (top right) and fO^** (bottom right) 
M(g yr~^, respectively. Each frame contains three curves corresponding, from left to right, to an atmosphere radius ratm = r^, 
ratm = lOrc and ratm = 5rc, respectively. When ratm = r^, the calculation is done for a protoplanet embedded in a standard 
disc model at a distance of 5 AU from the central star. The critical core mass is the mass beyond which no atmosphere at 
equilibrium can be supported at the surfa<;e of the core. 



